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There is little question that terrestrial ecosystems have a potentially pivotal role in 
determining both the direction and the rate of future changes in atmospheric CO; 
concentrations. There is some question, however, whether terrestrial ecosystems will 
contribute to global warming by releasing additional CO, into the atmosphere as a result 
of increasing respiration in a warming climate, or whether they will instead sequester 
additional carbon in response to the enhancing effects of atmospheric CO, on plant 
growth (Soloman and Cramer 1993). This latter response, the so-called CO, fertilization 
effect, has been advanced as an important negative feedback that has the potential for 
limiting increases in atmospheric CO» concentrations. The mechanism has often been 
invoked as an explanation of the 1.6-petagram (Pg = 1055 g)-carbon-per-year imbalance 
or "missing sink" in calculations of the global carbon budget. 

Despite the frequent incorporation of a CO, fertilization effect into ecosystem and 
global carbon models, questions concerning the magnitude of any additional storage of 
carbon in terrestrial biomass remain. On the one hand, direct experimental evidence 
shows that whole-plant carbon storage increases in response to elevated CO; concentra- 
tions (Kimball 1983); on the other, results from modeling exercises indicate that this 
enhanced short-term productivity of individual plants may not be sustained given limi- 
tations imposed on plant growth by the biogeochemical cycles of carbon and nitrogen 
(Rastetter et al. 1991). Historical tree ring chronologies are not definitive: in some 
instances tree growth has increased during the past century (LaMarche et al. 1984), 
whereas other chronologies indicate no such stimulation (Kienast and Luxmoore 1988). 
The issue is made more puzzling by the possibility that young or aggrading forests may 
be more likely than mature, closed-canopy forests to display increased productivity in 
response to elevated CO; concentrations (Graham et al. 1990). 
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There appears to be little hope that the uncertainty surrounding the magnitude of a 
CO, fertilization effect will be resolved soon. Progress is possible, however, using the 
results of controlled CO; exposure studies for species representative of ecosystems that 
store large amounts of carbon. Because forests are important reservoirs of carbon, a data 
base of over 390 prior observations on the response of trees to atmospheric CO; 
enrichment was compiled from 58 controlled-exposure studies. The data describe the 
growth of 73 boreal, temperate, and tropical tree species exposed to an approximate 
doubling of atmospheric CO;. This information was used to explore the magnitude of a 
CO, fertilization effect, to highlight species-specific differences in this response, and to 
evaluate whether a growth response to increasing CO; concentrations in forest trees 
might be constrained by limiting water and nutrients. Furthermore, we explored how the 
growth in controlled-exposure studies might be represented, applied, and interpreted 
within the context of global carbon models. Finally, we examined the potential conse- 
quences of a CO;-induced increase in plant productivity on litter decomposition and 
storage of carbon in soils. 


LITERATURE SEARCH, DATA PRESENTATION, AND ANALYSIS 


Literature describing the response of trees to CO; enrichment was identified through a 
computerized search of the AGRICOLA and BIOSIS data bases and by searching 
scientific journals. Information was also requested from investigators currently conduct- 
ing CO; enrichment studies. The data base is a comprehensive if not exhaustive repre- 
sentation of past and current investigations. 

Rather than simply reiterate the results for each of the studies, as has been the 
approach in previous reviews (Eamus and Jarvis 1989; Mousseau and Saugier 1992), 
we adopted instead the procedures of Kimball (1983), who, in documenting the 
growth response of agricultural crops to elevated CO» concentrations, looked not at 
the results of individual studies per se, but at responses characteristic of broad plant 
groups in a data base of over 430 observations. Although this type of analysis tends 
to ignore subtle differences between studies, it does offer the hope that the larger 
picture will be clearer and that generalities will appear that might otherwise remain 
obscure. 

~ The response of plant growth to increasing CO: concentration can be measured in 
many ways. À common technique is to calculate the relative change in total dry 
weight for plants grown at elevated CO; concentrations compared to that for their 
controls (Kimball 1983). In this approach a ratio of 1 would indicate no effect of 
elevated CO; on plant growth, whereas a ratio of 2 would indicate a doubling of dry 
mass in response to CO, enrichment. For the purposes of this presentation, relative 
growth responses were calculated for each observation. The data were transformed 
to logs to normalize the distribution and plotted as a frequency distribution showing 
the number of times that a given response was observed. The mean response for all 
observations, for the species-specific and biome-specific observations, and for the 
CO,-induced growth response of trees given limiting nutrient and water resources 
was calculated (Kimball 1983). In addition to the mean, 95% and 99.9% confidence 
intervals were also computed. 
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Figure 4.1. Frequency distribution of relative growth responses as calculated from the total 
dry mass of trees grown at ambient and elevated CO; concentrations. An average growth 
response for the entire data set was determined by averaging the log-transformed data to 
normalize the distribution. 


RESPONSE OF FOREST TREES TO CO; ENRICHMENT 


Based on the results of 58 controlled-exposure studies, the growth response of 73 tree 
species to elevated CO; concentrations was overwhelmingly positive (Figure 4.1). Of 
the 398 observations compiled in this study, 51 showed a relative growth response less 
than 1 (i.e., a value less than 0 on a log-transformed scale) and 31 showed values greater 
than 2 (i.e., logio > 0.30), indicating a twofold increase in plant dry mass with CO; 
enrichment. For the entire data set the relative growth response averaged 1.32, with a 
95% confidence interval of between 1.28 and 1.36 and a 99.9% confidence interval of 
between 1.25 and 1.39. A relative growth response of 1.32 is equivalent to a 32% 
increase in total plant dry mass (ie., above and below ground) in response to an 
approximate doubling of atmospheric CO» concentration. From a smaller number of 
Observations, this stimulation of whole-plant dry mass was found to be more or less 
equally partitioned between leaves (1.33), stems (1.29), and roots (1.38), but it was not 
a fair representation of canopy leaf area, which increased only marginally (1.13) in 
response to CO; enrichment. 

Separating the results by species indicated that the relative growth response varied 
fourfold across all 73 tree species, from below 1 to above 3 (Figure 4.2). The cause of 
the differences was not clear, except that species vary in response to CO; enrichment. 
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Figure 4.2. Species-specific estimates for the relative growth response as calculated from the 
total dry mass of trees grown at ambient and elevated CO; concentrations. Average growth 
responses were calculated by averaging the log-transformed data to normalize the distribu- 
tion. From left to right the species are: Populus trichocarpa, Quercus prinus, Betula sp., 
Myriocarpa longipes, Cecropia obtusifolia, Fraxinus lanceolata, Pseudotsuga menziesii, 
Populus euramericana, Carya ovata, Pinus virginiana, Senna multijuga, Liriodendron tulip- 
ifera, Platanus occidentalis, Picea sitchensis, Fraxinus americana, Acer saccharinum, Piper 
auritum, Nothofagus fusca, Pinus taeda, Quercus alba, Pinus strobus, Castanea sativa, Pinus 
radiata, Betula populifolia, Pinus ponderosa, Betula papyrifera, Populus trichocarpa x P. 
deltoides, Pinus echinata, Alnus rubra, Gliricidia sepium, Betula alleghaniensis, Picea 
mariana, Betula pendula, Tsuga canadensis, Pinus nigra, Pentaclethra macroloba, Prunus 
serotina, Robinia pseudoacacia, Pinus sylvestris, Acer rubrum, Trichospermum mexicanum, 
Pinus sp., Acer pensylvanicum, Salix x dasyclados, Picea glauca, Populus grandidentata, 
Acacia mangium, Pinus banksiana, Quercus rubra, Liquidambar styraciflua, Picea rubens, 
Alnus glutinosa, Picea pungens, Populus tremuloides, Elaeagnus angustifolia, Fagus syl- 
vatica, Cedrus atlantica, Populus deltoides, Juglans nigra, Fagus grandifolia, Picea sp., Abies 
fraseri, Acer saccharum, Picea abies, Ochroma lagopus, Pinus caribaea, Eucalyptus 
camaldulensis, Eucalyptus cypellocarpa, Tabebuia rosea, Eucalyptus pauciflora, Eucalyptus 
pulverulenta, Eucalyptus grandis, Pinus contorta. 


Those species most responsive to a doubling of atmospheric CO; were among the 
Eucalypts; estimates of the growth response for Eucalyptus grandis exceeded 3 and for 
the genus averaged over 2 (see the Appendix). No other genus-specific distinction in the 
growth response of forest trees to elevated CO, concentrations was apparent. Further 
separation of these species-specific results into those representative of three forest 
biomes indicated a trend: the relative growth response was greatest for trees of boreal 
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Table 4.1. Biome-Specific Estimates for the Relative Growth Response as 
Calculated from the Total Dry Mass of Trees Grown at Ambient and 
Elevated CO; Concentrations 


Lower limit Upper limit 
Forest biome -99.9% -95% Mean +95% +99.9% 
Boreal (n = 74) 1.24 1.30 1.38 1.46 1.53 
Temperate (n = 305) 1.24 1.27 1.31 1.35 1.39 
Tropical (n = 18) 0.98 1.10 1.25 1.41 1.54 


regions, intermediate for forests characterized by temperate zone broadleaf and conifer- 
ous species, and lowest for trees inhabiting tropical regions (Table 4.1). However, the 
differences were small and of questionable significance. 

Limitations imposed on the relative growth response by the availability of nitrogen 
and water were minimal (Table 4.2). Although the distinction between adequate and 
deficient as used here and in many of the experiments themselves is purely subjective 
(note the exception of Conroy et al. [1990a]), the growth response for forest trees 
exposed to CO; enrichment and grown under adequate supplies of nitrogen or water was 
1.32, a value only slightly higher than that for trees deficient in nitrogen or water 
(Table 4.2). In contrast, phosphorus deficiency did result in a lower, albeit still positive, 
relative growth response as compared with plants supplied with adequate phosphorus. 

The relative growth responses as estimated from the 58 controlled-exposure studies 
are representative of experiments in which CO, exposures lasted between 8 and 
798 days, with the average length of exposure being 136 days. The majority of these 
studies were conducted in growth chambers or greenhouses; a few were field studies 
using open-top chambers; most used transplanted seedlings grown in pots. Rarely were 
seedlings grown from seed under elevated CO; concentrations. Day and night tempera- 
tures ranged from 37/21°C to 20/10°C, with the 24-h average temperature being about 
22°C when pooled across all studies. No study considered the interactive effects of 
temperature on the growth response of trees to CO, enrichment. Ambient CO; concen- 
trations ranged from 300 to 400 ppm, with elevated concentrations of CO; ranging from 
415 to 2000 ppm (see the Appendix). With the possible exception of the work by Silvola 


Table 4.2. Estimates for the Relative Growth Response as Calculated from the Total Dry 
Mass of Trees Grown at Ambient and Elevated CO2 Concentrations Given Limiting 
Nutrient or Water Resources 


Lower limit Upper limit 

Resource —99.996 -95% Mean +95% +99.9% 
Nitrogen 

Adequate (n = 295) 1.24 1.27 1.32 1.36 1.40 

Deficient (n = 82) 1.16 1.22 1.30 1.38 1.45 
Phosphorus 

Adequate (n = 314) 1.25 1.28 1.33 1.37 1.40 

Deficient (n = 61) 1.09 1.14 1.21 1.29 1.35 
Water 

Adequate (n = 359) 1.26 1.28 1.32 1.37 1.40 


Deficient (n = 30) 1.16 1.21 1.28 1.35 1.41 
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and Ahlholm (1992), no study explicitly included a subambient CO; treatment. In the 
case of the Silvola and Ahlholm study, however, there was no ambient control. 


ESTIMATING A GLOBAL RESPONSE TO CO; ENRICHMENT 


Although the majority of controlled-exposure studies indicate that forest trees have the 
capacity to increase whole-plant carbon storage in response to enhanced concentrations 
of CO», many questions concerning the magnitude of such effects in natural ecosystems 
remain. Favoring the hierarchical development of models for describing ecosystem and 
global carbon dynamics, Reynolds et al. (1992) argue that the growth results obtained in 
controlled-exposure studies with plants must be interpreted with caution because they 
may not result in additional carbon storage in an ecosystem (see Kórner 1993). There 
are, of course, many processes that can intervene between the initial acquisition of 
carbon in photosynthesis, for which the response to rising CO; concentration is expected 
to be the greatest, and the ultimate expression of this in net ecosystem production. 
Integration of single-plant responses to heterogenous canopies and landscapes; inter- 
actions with limiting light, water, and nutrient resources; effects of plant life history and 
demography; competitive interactions among multiple species and life forms; and many 
more factors all need to be considered in extrapolating experimental results to the global 
carbon cycle. 

Although few would argue with the suggestion that intervening processes might 
indeed lessen the magnitude of a CO; fertilization effect as estimated from controlled- 
exposure studies (note the exception of Idso [1991]), it remains a difficult task to 
quantify the magnitude of a CO; fertilization effect at the ecosystem level, let alone on 
a global scale. Nevertheless, applying a CO) fertilization effect to ecosystem and global 
carbon models seems reasonable to mimic the assumed response of terrestrial vegetation 
to rising atmospheric CO; concentrations (Goudriaan and Ketner 1984; Kienast 1991; 
Gifford 1992; Polglase and Wang 1992). Here we describe, in some detail, how growth 
results obtained from controlled exposure studies can be cautiously used to derive 
estimates of the global CO; response by the terrestrial biosphere. 


Experimental Studies and a Global Biotic Growth Factor 


Keeling (1973) introduced a “biota growth factor," D, in his global carbon cycle model. 
The model represents carbon flux from the atmosphere to terrestrial biota and its 
response to increasing atmospheric CO; concentration. Keeling argued that “according 
to several studies in glass houses, the assimilation by plants at high CO» concentrations 
is roughly proportional to the logarithm of the ambient CO; concentration." In the form 
later used by Bacastow and Keeling (1973), this biota growth factor and its logarithmic 
dependency on CO; concentration were defined as 


Fav = Fool + BINN a/Nao)) (N/N 0) (4.1) 


where Fa» is the flux from the atmosphere to the terrestrial biota, F» is the initial 
preindustrial equilibrium flux from atmosphere to biota (and equivalently biota to 
atmosphere), N, and N, are the carbon contents of the atmosphere and terrestrial biota, 
respectively; and Nao and N; are their initial contents. The deficiencies of this log- 
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arithmic representation of CO; saturation have been frequently noted (e.g., Gifford 
1980; Gates 1985; Kohlmaier et al. 1987). Nevertheless, because of its historical prece- 
dence and its simplicity, Keeling's formulation has become one of the most common 
representations of terrestrial biospheric response to elevated CO; levels in global carbon 
cycle models. 

Keeling's biota growth factor applies to the global terrestrial biosphere as a whole, 
and its value has often, and appropriately, been determined by tuning for correspondence 
between observed and simulated changes in atmospheric CO; concentration. In its purest 
form, the global B of equation 4.1 cannot be determined experimentally. Nevertheless, 
the parameter D is explicitly related to plant response to CO; concentration, and it is 
reasonable to seek experimental evidence for appropriate a priori values of the biota 
growth factor. Within this context, some experimentalists have adopted B (often re- 
coined the "biotic growth factor") as a means of characterizing plant response to 
increasing CO, concentration (e.g., Gifford 1979a,b; Acock and Allen 1985; Allen et al. 
1987), and some global carbon cycle modelers have explicitly recognized an experimen- 
tally measured biotic growth factor that approximates the global parameter of the 
Keeling formulation (e.g., Kohlmaier et al. 1987). 

As noted in equation 4.1, Keeling's D applies to the global flux of carbon from 
atmosphere to biosphere (i.e, Fa). To estimate D from experimental studies of plant 
response to elevated CO; levels, we must translate between this global flux and experi- 
mental measurements. Net primary productivity (NPP) is normally used to approximate 
Fa and literature estimates of NPP are used to quantify NPP;. Thus we can rewrite 
equation (4.1) as 


NPP =NPP{1 + Bln(C./Ca0) (Cs /Coo) (4.2) 


where NPP, is the initial equilibrium NPP and C, and C, are carbon contents of the 
atmosphere and biosphere, respectively. From equation 4.2 we can define D as 


_ (NPP — NPP») | NPPo 


ps In (Ca / Cao) i 


If we assume that NPP is independent of the mass of the terrestrial biosphère (e.g., 
Goudriaan and Ketner 1984; Kohlmaier et al. 1987) then 


NPP = NPP,[1 + Biln(C,/C.0)] (4.4) 
and 


_ (NPP —NPPo»)/ NPPo 


B. In (C, / Cao) e 


For small changes in the terrestrial biosphere By = Bı, and for a monotonically 
increasing terrestrial biosphere Bo < By. 
For small changes in atmospheric CO; concentration, 


B- (NPP — NPP») / NPPo 
UU Gi= Ca)! Ca 


We identify this approximation as B». Thus for small changes in the terrestrial biosphere 
and atmospheric CO», B. = B; = Bo, and D is approximately the fractional increase in NPP 


(4.6) 
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for a unit increase in atmospheric CO, concentration. B; of equation 4.6 is often 

described as the biotic growth factor and ascribed to Bacastow and Keeling (1973). We 

note that it is in fact an approximation of the B (equation 4.3) described in that paper. 
Net primary productivity can be defined as 


NPP =AB+L+C (4.7) 


where AB is the change in vegetation biomass over some period of time and L and C 
are litter production and the loss of vegetation biomass to herbivore consumption, 
respectively, over the same period of time. If we assume that herbivory is absent or 
negligible in experimental CO, exposure studies, and that any litter produced is 
accounted for as a change in biomass, then we can define P as AB over the period of 
the experimental exposure and equate NPP with P. Then, if we assume that the 
functional relationship of equation 4.2 among global NPP, biospheric carbon, and 
atmospheric CO; changing over tens to hundreds of years also characterizes, at least 
approximately, the relationship among plant biomass production, plant biomass, and 
relatively short-term exposure to elevated CO; levels under controlled experimental 
conditions, we can define the relationship 


P.=P.{1 + BIn(C./C.)KB./B.) (4.8) 


where P. is plant biomass production at control CO; concentrations C., P. is production 
at experimental or elevated CO, concentrations C., and B. and B, are plant biomass at 
control and elevated CO, concentrations, respectively. 

From equation 4.8 and the corresponding translations of equations 4.5 and 4.6, we 
can define the following experimental growth factors: 


_ IP. (G3. B.) - P. ]/ P. 


BE; In (C, / C.) (4.92) 
_(@.-P.)1 Pe 

BE, = In (C./ C) (4.9b) 
_(@.—P.)/ P. 

E mA ud 


Upon equating P with NPP, BE; of equation 4.9c is equivalent to the experimental 
fertilization factor B of Kohlmaier et al. (1987). 

Controlled CO, exposure studies often report only the biomass of treatments at the 
end of the experiment and not the explicit biomass production P, analogous to NPP, 
called for by equation 4.9. We have defined P as AB over the period of the CO; exposure 
experiment, or 


P. =Ba -Bo (4.10a) 
P.=Ba-Ba (4.10b) 


where t is the length of exposure and the subscript 0 indicates initial values at t= 0 such 
that Ba, for example, is the biomass of the control exposure at the end of the experiment. 
B. and B. of equation 4.8 are defined as B« and B,, respectively. 

If the experimental treatment begins with seed germination, as might be the case with 
annual crop plants, the Bos are essentially zero, P = B,. BE of equation 4.9a is ill defined 
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as zero by this substitution. However, BE, and BE of equations 4.9b and 4.9c are 
unaffected by substituting E, for P. If the initial biomass of the experimental material is 
nonzero, but appreciably the same for each treatment, as is the case for most studies with 
woody perennials, then the experimental growth factors can be evaluated as 


ap. Ba / Be)(Ber — Ba) — (B4 — Boo)| / Ba — Bo 


u= In (C./ C) MAI 
_ (Ba - Bo) / (Ba - Ba) 
OMNE TOIT M (4.11b) 


ap _ Ba- Bo)! (Ba — Boo) 


Mae -CC i 


To the extent that Bo # B.o, the BEs of equations 4.11 are approximations of the BE of 
equations 4.9. If B.o > Ba they are overapproximations; they are underapproximations if 
Ba < Bo. 

The relative growth response used earlier in this chapter to summarize the controlled 
exposure studies was defined as 


P = BalBe (4.12) 


where, for the entire data set, p = 1.32. For Ba much greater than Bao, we can make the 
following approximations: 


ac [p B. Ba) Bal 1 

OMNE T RT (4.132) 
E acp - 
Mime (4.13b) 
ĝe = Pa (4.130) 


and estimate a corresponding biotic growth factor for each approximation 
(Table 4.3). Since B.o = Bao < Ba € Ba BEo is again poorly approximated as being 
close to zero. 


Table 4.3. Estimates of a Global Biotic Growth Factor Based 
on the Results of Controlled-Exposure Studies 


i NBEja [20 BEic Bid 
I 70.46 0.46€ 0.46 0.46ye € 0.46 
2 =0.32 0.32e 0.32€ 0.32ye < 0.32 


"Estimates of N BE based on equations 4.14. 

Estimates of ÎE of equations 4.11 based on the approximation of equations 4.14, where 
0 «€ € 1, and probably close to 1. . 

“Estimates of BE of equations 4.9 from estimates of BE, assuming initial biomass of 
control and elevated CO» treatments is equal. 

"Estimates of the global B of equations 4.3, 4.5, and 4.6 by equations 4.15, where 0 < 
ysl. 
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The control and elevated CO; treatments vary among the studies summarized in the 
Appendix, but we can normalize for a doubling of CO; concentration (C, = 2C.) such 
that 


Ap [p (Ba Ba) / Bal - 1 
NBEo = "r5 (4.142) 
ae p-1 
Nh Po (4.14b) 
NBE) = p—1 (4.14c) 


Once again, NBEo is poorly approximated as close to zero. When the elevated CO; 
treatment is greater than twice the control, these normalized approximations NBE, and 
NBE: are overestimates of the BE. We also note that for a doubling of atmospheric CO; 
concentration, the B; approximation of B, underestimates B; by roughly 30%. 

Without data on rates of production (P) or the initial biomass (Bp) of control and 
elevated treatments (data absent from most reports of controlled CO; exposure studies, 
including those summarized here), the experimental growth factor BE, (or BEo) cannot 
be appropriately evaluated. However, estimates of the 2xCO; normalized NBE, and 
NBE, for the mean relative growth response of the entire data set presented here (p = 
1.32) are given in Table 4.3. As noted previously, these estimates are approximations of 
the BE of equations 4.11 and are probably overestimates. The correction for this approx- 
imation is indicated by £ (Table 4.3), where 0 < £ € 1. Because we believe B, to be much 
less than B« in most of the studies summarized here, and most elevated CO; concentra- 
tions are at least approximately a doubling of the control concentrations, we believe that 
g is close to 1. We assume that Bo = B4; and accept the approximations of BE as 
estimates of the BE of equation 4.9. 

In translating from the experimental growth factors to an estimate of the global B, we 
make the following definitions: 


Bo = YBEs (4.15a) 
Bi. = YBE; (4.15b) 
Bo = YBE2 (4.15c) 


where y is a proportionality coefficient that incorporates the multitude of assump- 
tions, errors, and uncertainties involved in translating from the measurements of 
short-term experimental exposures of small (generally young) plants in controlled 
chambers or small field exposures to a model parameter of long-term CO; response 
by the terrestrial biosphere. The value of y is highly uncertain. Experimental evi- 
dence is mixed, but there is some indication that as the level of hierarchical system 
organization increases (e.g., from leaf photosynthetic rates to whole-plant growth), 
the CO; response declines (e.g., Allen et al. 1987) and that with increased exposure 
time acclimation can occur with a decline in the magnitude of the CO; response 
(Oechel and Strain 1985). Additional experiments, modeling, analysis, and theoreti- 
cal treatment are required to substantiate or refute these trends, but for now we 
assume that 0 € y < 1. The resulting estimates of global B as defined in equations 4.3, 
4.5, and 4.6 are given in Table 4.3. 
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Carbon Sequestration and Soil Organic Matter 


The capacity of forest trees to increase whole-plant carbon storage is only one of many 
components that will dictate whether similar increases in carbon sequestration will be 
observed in natural ecosystems, where carbon will be stored not only in biomass but also in 
soils. In many terrestrial ecosystems, for example, the amount of carbon sequestered in soil 
organic matter greatly exceeds that in living plants. Similarly, on the global scale over twice 
as much carbon is stored within soil organic matter and its component fractions as in plant 
biomass (Post et al. 1990). Given these distinctions, if we assume for the moment that plant 
growth will increase in response to rising CO, concentrations, then the emphasis shifts away 
from one that focuses strictly on vegetation toward one of understanding how this increased 
biomass, once shed from the plant, decays and becomes soil organic matter. These are the 
processes that will largely govern the rate and magnitude of long-term carbon storage under 
changing atmospheric CO) concentrations. 

The amount of carbon stored in soils is determined by the rate of organic matter input 
through litterfall (i.e., leaves, branches, boles), root exudates, fine root turnover, and the 
rate of decay. However, a large fraction of soil organic matter has a long turnover time 
(around 20-50 years), and this introduces a considerable time response for sequestering 
in soils any additional carbon that may result from increasing plant productivity owing 
to CO, enrichment. If, for example, the rate of organic matter input into soil increases by 
30% in response to a doubling of atmospheric CO, concentration, then, all other factors 
remaining the same, soil carbon should also increase by an equivalent amount. An 
important question to ask in this regard, however, is how long it will take for this new 
equilibrium to become established. 

Using the Rothamsted turnover model (Jenkinson 1990), it is possible to explore the 
time lags involved in soil carbon sequestration and to link CO;-induced effects on plant 
productivity with decomposition. Plant residues are divided into two components. The 
first includes simple compounds subject to rapid uptake, transformation, and mineraliza- 
tion by decomposers. This pool is decomposable plant matter (DPM) and has a turnover 
time of 0.1 year. The second pool contains components such as lignin that are not readily 
attacked; it is identified as recalcitrant plant matter (RPM), with a tumover time of 
3 years. Soil humus is divided into three components with turnover times differing by an 
order of magnitude. The soil biomass (BIO, turnover time 1.5 years) component is 
mostly live microbes and microbial products such as extracellular enzymes, whereas the 
active soil humus (HUM, turnover time 50 years) component consists of organic com- 
pounds that are protected physically or that exist in chemical forms that are biologically 
resistant to decomposition. Finally, the inert soil organic matter (IOM) component 
consists of chemically refractory and physically protected forms with turnover times in 
excess of 1000 years. In the simulations presented here, the IOM is truly inert, neither 
increasing nor decreasing during the simulations. Following Jenkinson et al. (1991) we 
set this pool at 3 tons/hectare. 

The fate of soil carbon as estimated from the Rothamsted model was investigated for 
a hypothetical 200-year period during which organic matter inputs to the soil were 
gradually increased by 3096 (Figure 4.3). Modeled results using Rothamsted Experiment 
Station climatic conditions (mean annual temperature 9.3*C) showed that the litter pools 
(DPM and RPM) and decomposer biomass (BIO) responded quickly, because of their 
rapid turnover times, and tracked closely the projected increases in organic matter inputs 
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Figure 4.3. Soil carbon responses to increasing inputs as simulated by the Rothamsted 
turnover model. Two scenarios are presented, one using the climatic conditions of 
Rothamsted Experiment Station (A), where annual temperatures average 9.3°C, and another 
using the warmer conditions of Manaus, Brazil (B), where annual temperatures average 
26.9°C. Organic matter inputs were assumed to increase linearly by 30% over a 200-year 
period beginning in 1860. The DPM pool increased at the same rate as total inputs. The 
RPM and BIO pools lagged closely behind the organic matter increases, although RPM in 
the Manaus simulations did equilibrate as readily as total inputs and overlay the DPM line 
in this subplot. The HUM pool lagged behind the other organic matter components and 
resulted in the response of the litter and active soil pools (DPM + RPM + BIO + HUM) not 
reaching equilibrium during the simulation period. The magnitude of the lag was related to 
the rate of decomposition, which was considerably higher with Manaus climate than with 
Rothamsted climate. The IOM pool of the Rothamsted model did not change during the 
simulation, is not plotted here, and was not used in the calculation of the total litter and soil 
response. 
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(Figure 4.3A). Increased inputs because of enhanced plant productivity were imme- 
diately reflected in these component fractions. These pools, however, contained less than 
30% of total soil organic matter. The much larger active soil humus (HUM) pool, which 
comprises 70% of the soil organic matter, increased at a much slower rate. As a result, 
the total increase in soil organic matter after 200 years was only 17%, substantially less 
than one might expect from a 30% increase in litter inputs. Nonetheless, such an increase 
resulted in incorporating 15% of the enhanced productivity over this 200-year period 
into soil organic matter, mostly as active humus. Model runs using the warmer climatic 
conditions for Manaus, Brazil (mean annual temperature 26.9°C) resulted in higher 
decomposition rates (Figure 4.3B). The HUM pool, which in this case comprised 34% of 
the total soil carbon, increased faster under these conditions than for the Rothamsted 
scenario. As a result, the total soil carbon increased by 25% over the 200-year period, 
although much less of the enhanced biomass production during this time was seques- 
tered, with only 4% retained in soil carbon pools—30% in RPM and 67% in HUM. 
From these modeled results we observe that the storage of carbon in biomass and in 
soils may not be tightly coupled: an increase in whole-plant productivity resulting from 
rising CO, concentrations may not necessarily lead to a proportional increase in soil 
carbon storage, despite an equivalent increase in litterfall. There are two reasons for this. 
First, the turnover times of soil organic matter and its component fractions are suffi- 
ciently long that any increased input in the form of litterfall and its subsequent incorpo- 
ration into soil organic matter will not occur immediately, but will be delayed depending 
on the turnover times of individual soil carbon pools and decomposition rates. Second, 
as indicated by the contrasting results of the Rothamsted and Manaus climatic effects, 
increasing decomposition rates will increase the speed at which the slower organic 
matter pools approach equilibrium, but as temperatures increase less carbon will be 
sequestered in soil carbon pools and will be returned to the atmosphere as CO2. 


GLOBAL CARBON MODELS AND A CO; FERTILIZATION EFFECT 


As alluded to earlier, global carbon models frequently mimic the CO, fertilization effect in 
terrestrial vegetation by incorporating a biotic growth factor. In an interesting application of 
the biotic growth factor, Polglase and Wang (1992) divided the globe into 10 representative 
biomes and then used the Rothamsted soil carbon model (Jenkinson 1990)—along with 
estimates of areal distribution; net primary productivity; allocation of biomass to leaves, 
branches, stems, and roots (each with a characteristic residence time); and a reference-state 
biotic growth factor of 0.3—to calculate biome-specific and global estimates of CO,- 
enhanced carbon storage. The biotic growth factor was made a function of annual air 
temperature, with the growth factor being lower in cooler and higher in warmer climates. 
They applied this growth factor to an NPP function like that of equation 4.4. Their results 
indicated that the total CO;-enhanced carbon storage in 1990 was 1.14 Pg/year. The most 
important biomes for sequestering carbon were tropical forests, savannas, and temperate 
zone forests. Of the globally sequestered carbon, 42% was stored in vegetation, 19% in litter, 
and 39% in soil organic matter. Our estimate of a larger biotic growth factor, B, € 0.46, would 
presumably result in even greater CO;-enhanced carbon storage. 

Calculating probable rates of global carbon sequestration caused by land-use change, 
primarily deforestation in the tropics, and a CO. fertilization effect, Taylor and Lloyd (1992) 
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Table 4.4. Summary of Studies That Have Incorporated a Biotic Growth 
Factor into Models of Global Carbon Sequestration 


Biome Biotic growth Rate of carbon 

Reference divisions factor storage (Pg/year) 
Oeschger et al. (1975) 1 0.20 — 
Gifford (1980) 9 0.60 1.0 
Goudriaan and Ketner (1984) 6 0.50 d 
Kohlmaier et al. (1987) 4 0.25-0.50 14 
Kohlmaier et al. (1989) — 0.38 -— 
Gifford (1993) 1 0.36 23 
Polglase and Wang (1992) 10 0.30 14 

10 0.02-0.54 L5 

10 0.02-0.54 13 
Taylor and Lloyd (1992) 14 0.36 24 


Note: Because of the many assumptions used in each model the rate of carbon storage cannot be ascribed 
solely to a CO» fertilization effect. 


similarly invoked a biotic growth factor of 0.37 (i.e., a 2596 increase in NPP associated 
with a doubling of atmospheric CO, concentration) for 14 major ecosystem types. They 
suggested that the total net sink effect in 1985 for carbon was 2.4 Pg/year globally, with 1996 
of the terrestrial carbon sink in the boreal zone, 37% in the temperate zone, and 44% in the 
tropics. Their biotic growth factor corresponds to our f» (i.e., 0.32) shown in Table 4.3. Our 
estimate of this biotic growth factor is only slightly smaller than that used by Taylor and 
Lloyd and would presumably reduce the estimates of a net sink only slightly. 

The examples cited, along with those obtained by Gifford (1993) using the C-Quest 
model, and a sampling of other studies that include a biotic growth factor in global 
carbon models (Table 4.4) indicate that through a CO; fertilization effect on net primary 
productivity, terrestrial ecosystems can sequester vast amounts of carbon and that this 
process might indeed function as an important negative feedback on future increases in 
global CO; concentrations. For example, in simulating the concentrations of atmo- 
spheric CO, for the year 2030, Goudriaan and Ketner (1984) indicated that by reducing 
the biotic growth factor from 0.50 to 0.25, and thereby lessening the negative feedback 
potential of terrestrial vegetation, atmospheric CO; concentrations would increase by an 
estimated 21 ppm. Likewise, Polglase and Wang (1992) also point out the sensitivity of 
carbon sequestration to the biotic growth factor, suggesting that a twofold increase in B 
gives rise to a doubling of ecosystem and global CO;-enhanced carbon storage. Given 
this sensitivity of model estimates to the often subjective choice of a biotic growth 
factor, it is important to consider carefully and to understand fully the assumptions and 
translations involved in using controlled-exposure studies to estimate a global biotic 
growth factor, and how that estimate is applied in global carbon models. 


CONCLUSIONS 


Where do the results obtained from controlled-exposure studies begin to provide 
answers for questions of truly global proportions? So far the link between results 
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obtained fromexperimental studies and the application of suchresults on the ecosystem 
or global scale has been a limited one. This has not, however, been the case because 
of a lack of effort, Rather, it has occurred because of the realization that answers to 
ecosystem and global concerns may not always lend themselves to direct experimen- 
tation, and as a result modelers have been left, in many cases, with little more than 
intuition to guide them. This comment is not intended to be derogatory; instead it 
suggests that experimental studies, among all the important questions that they are 
addressing, may not be providing the simple information needed to verify or refute 
even the most basic of assumptions required in global carbon models. Within the 
context of the CO, fertilization effect and the data summarized herein, four points 
are illustrative: 


1. Despite the fact that boreal, temperate, and tropical forests occupy regions of the 
globe with vastly different temperature regimes, these differences are not reflected 
in the growth conditions chosen for CO; enrichment studies. Of the 58 controlled- 
exposure studies considered here, virtually all were conducted at temperatures of 
20—25'C, temperatures that apparently were chosen with little thought to the species” 
natural range. 

2. Even if one wanted to scale the CO; fertilization effect for changing temperatures, 
we know very little about what such a correction curve would look like. Although 
Polglase and Wang (1992) derive a solution whereby B can be corrected for 
temperature, there has to date been no attempt to examine the potential for a 
temperature x CO; level interaction in forest trees. Because of the presumed 
association between rising CO, concentrations and global warming, there is an 
obvious need to pursue this area of research. 

3. It is clear that to predict future atmospheric concentrations of CO; and hence to 
evaluate the potential for global warming, modelers must first match the predic- 
tive capabilities of a carbon model against the historical record. The period 
considered in this chapter was characterized by CO; concentrations below present- 
day levels, but no study has examined the growth response of long-lived woody 
perennials to subambient CO; concentrations. 

4. Given the need to equate experimentally derived estimates of biomass or 
biomass productivity with net primary productivity, it is unfortunate that so 
few studies report an initial plant dry mass at the beginning of the CO; 
exposure. Doing so would help eliminate some of the assumptions, approxima- 
tions, and uncertainties in estimating a global B from controlled-exposure 
studies. 


Clearly, experimental studies have done much to address the needs of modeling the 
global carbon cycle, yet it is equally clear that much work remains ahead. In assessing 
the capacity of terrestrial vegetation to serve as an important biospheric feedback to 
rising CO; concentrations and associated climate change, it is apparent that we have only 
begun to meet the challenges. By working together, experimentalists and modelers must 
strive to complement each others' strengths and weaknesses to define and interpret 
empirical and mechanistic descriptions of terrestrial biota and subsequently incorporate 
them into ecosystem and global carbon models. 
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APPENDIX: SPECIES-SPECIFIC ESTIMATES FOR THE RELATIVE GROWTH 
RESPONSE AS CALCULATED FROM THE TOTAL DRY MASS OF TREES 
GROWN AT AMBIENT AND ELEVATED CO; CONCENTRATIONS 


Ambient and elevated concentrations of CO; used in each study varied, as did the 
duration of each study, and this information is provided for the reader's convenience. 
NA denotes data not available. 


Growth CO2 Duration of Relative 


concentration experiment growth 
Species (ppm) (days) response Reference 
Abies fraseri 362,713 365 1.76 Samuelson and Seiler (1992) 
Acacia mangium 354,711 95 1.40 Ziska et al. (1991) 
Acer pensylvanicum 350, 700 165 1.36 Bazzaz and Miao (1993) 
Acer rubrum 400, 700 60 1.05 Bazzaz et al. (1990) 
350, 700 70 1.58 Bunce (1992) 
350, 700 66 1.36 Miao et al. (1992) 
350, 700 165 1.34 Bazzaz and Miao (1993) 
Acer saccharinum 300, 600, 1200 35 1.75 Carlson and Bazzaz (1980) 
350, 500, 700 90 0.87 Williams et al. (1986) 
350, 700 55 1.66 Bunce (1992) 
Acer saccharum 400, 700 100 2.52 Bazzaz et al. (1990) 
400, 800 85 1.55 Nobel et al. (1992) 
Alnus glutinosa 350, 700 98 1.49 Norby (1987) 
Alnus rubra 350, 650 1547 1.27 Amone and Gordon (1990) 
Betula alleghaniensis 350, 700 165 1.27 Bazzaz and Miao (1993) 
Betula papyrifera 400, 700 60 1.23 Bazzaz et al. (1990) 
Betula pendula 350, 700 28-42 1.31 Pettersson and McDonald 
(1992) 
Betula populifolia 350, 700 66 1.30 Miao et al. (1992) 
350, 700 165 1.17 Bazzaz and Miao (1993)) 
Betula sp. 350, 2000 40 0.93 Hardh (1966) 
Carya ovata 350, 500, 700 90 1.08 Williams et al. (1986) 
Castanea sativa 350, 700 365 1.23 El Kohen et al. (1992) 
350, 700 NA 1.14 Mousseau (1992) 
350, 700 730 1.14 Mousseau and Saugier (1992) 
350, 700 365 1.20 El Kohen et al. (1993) 
Cecropia obtusifolia 350, 525, 700 11 1.00 Reekie and Bazzaz (1989) 
Cedrus atlantica 350, 800 274 1.66 Kaushal et al. (1989) 
Elaeagnus angustifolia 350, 700 84 1.61 Norby (1987) 
Eucalyptus 330, 660 84 221 Wong et al. (1992) 
camaldulensis 
Eucalyptus 330, 660 84 2.48 Wong et al. (1992) 
cypellocarpa 
Eucalyptus grandis 340, 660 42 3.34 Conroy et al. (1992) 
Eucalyptus pauciflora 330, 600 100 2.75 Wong et al. (1992) 
Eucalyptus 330, 660 98 2.88 Wong et al. (1992) 
pulverulenta 
Fagus grandifolia 400, 700 60 1.74 Bazzaz et al. (1990) 
Fagus sylvatica 350, 700 365 1.62 El Kohen et al. (1993) 
Fraxinus americana 350, 700 165 1.19 Bazzaz and Miao (1993) 
Fraxinus lanceolata 350, 500, 700 90 1.03 Williams et al. (1986) 


Gliricidia sepium 350, 650 31-71 1.26 Thomas et al. (1991) 
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Growth CO; Duration of Relative 


concentration experiment growth 
Species (ppm) (days) response Reference 
Juglans nigra 325, 1500 NA 1.46 Tinus (1976) 
Liquidambar 340,910 91 1.60 Rogers et al. (1983) 
styraciflua 
350, 675, 1000 8-113 1.48 Tolley and Strain (19842) 
350, 675, 1000 54-78 1.35 Tolley and Strain (1984b) 
350, 500, 700 224 1.44 Sionit et al. (1985) 
Liriodendron tulipifera 350, 500, 700 90 0.95 Williams et al. (1986) 
367, 692 168 1.72 O'Neill et al. (1987) 
371, 493, 787 168 1.15 Norby and O'Neill (1991) 
Myriocarpa longipes 350, 525, 770 111 0.96 Reekie and Bazzaz (1989) 
Nothofagus fusca 340, 640 120 1.17 Hollinger (1987) 
Ochroma lagopus 350, 675 60 1.79 Oberbauer et al. (1985) 
Pentaclethra 350, 675 123 1.30 Oberbauer et al. (1985) 
macroloba 
Picea abies 350, 900 21 1.78 Yeatman (1970) 
Picea glauca 350, 900 21 1.61 Yeatman (1970) 
350, 700 30-100 1.37 Brown and Higginbotham 
(1986) 
Picea mariana 340, 1000 21-42 1.27 Campagna and Margolis 
(1989) 
Picea pungens 325, 1200 243-365 1.50 Tinus (1972) 
Picea rubens 350, 700 122 1.12 Shipley et al. (1992) 
374, 713 365 1.66 Samuelson and Seiler (1994) 
362, 711 152 1.49 Samuelson and Seiler (1993) 
Picea sitchensis 350, 700 213-426 1.14 Lee et al. (1992) 
Picea sp. 350, 2000 40 1.75 Hardh (1966) 
Pinus banksiana 350, 900 21 1.40 Yeatman (1970) 
Pinus caribaea 340, 600 343 2.07 Conroy et al. (1990b) 
Pinus contorta 300, 1000 152 5.20 Higginbotham et al. (1985) 
Pinus echinata 368, 695 238-287 145 Norby et al. (1987) 
360, 700 42-168 1.36 O'Neill et al. (1987b) 
Pinus nigra 350, 800 274 1.29 Kaushal et al. (1989) 
Pinus ponderosa 325, 1200 243-365 1.35 Tinus (1972) 
355, 530, 705 122-243 1.28 Ball et al. (1992) 
Pinus radiata 330, 660 63-154 131 Conroy et al. (19862) 
330, 660 63-154 1.20 Conroy et al. (1986b) 
340, 640 120 1.27 Hollinger (1987) 
330, 660 154 1.24 Conroy et al. (1988) 
340, 660 798 1.19 Conroy et al. (1990a) 
Pinus strobus 400, 700 100 1.19 Bazzaz et al. (1990) 
Pinus sylvestris 350, 900 21 1.32 Yeatman (1970) 
Pinus taeda 350, 675, 1000 8-84 1.15 Tolley and Strain (1984a) 
350, 675, 1000 54-82 1.08 Tolley and Strain (1984b) 
350, 500, 700 224 1.58 Sionit et al. (1985) 
350, 500, 650 155-512 1.24 Ball et al. (1992) 
350, 700 181 1.42 T. J. Tschaplinski et al. (1993) 
Pinus virginiana 340, 415, 640, 122 1.09 CDIC (1985) 
940 
Pinus sp. 350, 2000 40 1.34 Hardh (1966) 
Piper auritum 350, 525, 700 111 1.16 Reekie and Bazzaz (1989) 


(Table continued on next page) 
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Growth CO; Duration of Relative 
concentration experiment growth 
Species (ppm) (days) response Reference 
Platanus occidentalis 300, 600, 1200 35 1.22 Carlson and Bazzaz (1980) 
350, 500, 700 90 1.09 Williams et al. (1986) 
Populus deltoides 300, 600, 1200 35 1.70 Carlson and Bazzaz (1980) 
Populus euramericana 350, 700 92 1.72 Radoglou and Jarvis (1990) 
350, 700 228 0.65 R. Ceulemans (unpublished 
data) 
Populus grandidentata 361, 707 70 1.27 Curtis and Teeri (1992) 
342,692 152 1.50 Zak et al. (1993) 
Populus tremuloides 350, 750 30-100 1.68 Brown and Higginbotham 
(1986) 
Populus trichocarpa 350, 700 92 1.12 Radoglou and Jarvis (1990) 
350, 700 228 0.64 R. Ceulemans (unpublished 
data) 
Populus trichocarpa 350, 700 92 1.24 Radoglou and Jarvis (1990) 
x P. deltoides 
350, 700 228 1.23 R. Ceulemans (unpublished 
data) 
Prunus serotina 400, 700 60 1.32 Bazzaz et al. (1990) 
Pseudotsuga menziesii 340, 640 120 1.03 Hollinger (1987) 
Quercus alba 360, 700 42-210 1.62 O'Neill et al. (1987a) 
386, 496, 793 35-268 1.07 Norby and O"Neill (1989) 
Quercus rubra 350, 500, 700 90 1.09 Williams et al. (1986) 
350, 700 165 1.98 Bazzaz and Miao (1993) 
Quercus prinus 350, 700 64 0.89 Bunce (1992) 
Robinia pseudoacacia 350, 700 105 1.32 Norby (1987) 
Salix x dasyclados 300, 500, 700 122 1.36 Silvola and Ahlholm (1992) 
Senna multijuga 350, 525, 700 111 1.10 Reekie and Bazzaz (1986) 
Tabebuia rosea 354,711 95 2.63 Ziska et al. (1991) 
Trichospermum 350, 525, 700 111 1:35. Reekie and Bazzaz (1989) 
mexicanum 
Tsuga canadensis 400, 700 100 1.28 Bazzaz et al. (1990) 
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